Frequency stabilization and noise-induced spectral narrowing in resonators with zero dispersion.

Mechanical resonators are widely used as precision clocks and sensitive detectors that rely on the stability of their eigenfrequencies. The phase noise is determined by different factors including thermal noise, frequency noise of the resonator and noise in the feedback circuitry. Increasing the vibration amplitude can mitigate some of these effects but the improvements are limited by nonlinearities that are particularly strong for miniaturized micro- and nano-mechanical systems. Here we design a micromechanical resonator with non-monotonic dependence of the eigenfrequency on energy. Near the extremum, where the dispersion of the eigenfrequency is zero, the system regains certain characteristics of a linear resonator, albeit at large amplitudes. The spectral peak undergoes narrowing when the noise intensity is increased. With the resonator serving as the frequency-selecting element in a feedback loop, the phase noise at the extremum amplitude is ~3 times smaller than the minimal noise in the conventional nonlinear regime.

Time-resolved measurement of Landau-Zener tunneling in periodic potentials.

We report time-resolved measurements of Landau-Zener tunneling of Bose-Einstein condensates in accelerated optical lattices, clearly resolving the steplike time dependence of the band populations. Using different experimental protocols we were able to measure the tunneling probability both in the adiabatic and in the diabatic bases of the system. We also experimentally determine the contribution of the momentum width of the Bose condensates to the temporal width of the tunneling steps and discuss the implications for measuring the jump time in the Landau-Zener problem.

Maximal width of the separatrix chaotic layer.

The main goal of the paper is to find the absolute maximum of the width of the separatrix chaotic layer as function of the frequency of the time-periodic perturbation of a one-dimensional Hamiltonian system possessing a separatrix, which is one of the major unsolved problems in the theory of separatrix chaos. For a given small amplitude of the perturbation, the width is shown to possess sharp peaks in the range from logarithmically small to moderate frequencies. These peaks are universal, being the consequence of the involvement of the nonlinear resonance dynamics into the separatrix chaotic motion. Developing further the approach introduced in the recent paper by Soskin [Phys. Rev. E 77, 036221 (2008)], we derive leading-order asymptotic expressions for the shape of the low-frequency peaks. The maxima of the peaks, including in particular the absolute maximum of the width, are proportional to the perturbation amplitude times either a logarithmically large factor or a numerical, still typically large, factor, depending on the type of system. Thus, our theory predicts that the maximal width of the chaotic layer may be much larger than that predicted by former theories. The theory is verified in simulations. An application to the facilitation of global chaos onset is discussed.

Matching of separatrix map and resonant dynamics, with application to global chaos onset between separatrices.

We have developed a general method for the description of separatrix chaos, based on the analysis of the separatrix map dynamics. Matching it with the resonant Hamiltonian analysis, we show that, for a given amplitude of perturbation, the maximum width of the chaotic layer in energy may be much larger than it was assumed before. We use the above method to explain the drastic facilitation of global chaos onset in time-periodically perturbed Hamiltonian systems possessing two or more separatrices, previously discovered [S. M. Soskin, O. M. Yevtushenko, and R. Mannella, Phys. Rev. Lett. 90, 174101 (2003)]. The theory well agrees with simulations. We also discuss generalizations and applications. The method may be generalized for single-separatrix cases. The facilitation of global chaos onset may be relevant to a variety of systems, e.g., optical lattices, magnetic and semiconductor superlattices, meandering flows in the ocean, and spinning pendulums. Apart from dynamical transport, it may facilitate noise-induced transitions and the stochastic web formation.

Langevin stabilization of molecular-dynamics simulations of polymers by means of quasisymplectic algorithms.

Algorithms for the numerical integration of Langevin equations are compared in detail from the point of view of their accuracy, numerical efficiency, and stability to assess them as potential candidates for molecular-dynamics simulations of polymeric systems. Some algorithms are symplectic in the deterministic frictionless limit and prove to stabilize long time-step integrators. They are tested against other popular algorithms. The optimal algorithm depends on the main goal: accuracy or efficiency. The former depends on the observable of interest. A recently developed quasisymplectic algorithm with great accuracy in the position evaluation exhibits better overall accuracy and stability than the other ones. On the other hand, the well-known BrunGer-Brooks-Karplus [Chem. Phys. Lett. 105, 495 (1982)] algorithm is found to be faster with limited accuracy loss but less stable. It is also found that using higher-order algorithms does not necessarily improve the accuracy. Moreover, they usually require more force evaluations per single step, thus leading to poorer performances.

Divergence of the chaotic layer width and strong acceleration of the spatial chaotic transport in periodic systems driven by an adiabatic ac force.

We show for the first time that a weak perturbation in a Hamiltonian system may lead to an arbitrarily wide chaotic layer and fast chaotic transport. This generic effect occurs in any spatially periodic Hamiltonian system subject to a sufficiently slow ac force. We explain it and develop an explicit theory for the layer width, verified in simulations. Chaotic spatial transport as well as applications to the diffusion of particles on surfaces, threshold devices, and others are discussed.

Solution of the boundary value problem for optimal escape in continuous stochastic systems and maps.

Topologies of invariant manifolds and optimal trajectories are investigated in stochastic continuous systems and maps. A topological method is introduced that simplifies the solution of boundary value problems: The activation energy is calculated as a function of a set of parameters characterizing the initial conditions of the escape path. The method is applied explicitly to compute the optimal escape path and the activation energy for a variety of dynamical systems and maps.

Quasisymplectic integrators for stochastic differential equations.

Two specialized algorithms for the numerical integration of the equations of motion of a Brownian walker obeying detailed balance are introduced. The algorithms become symplectic in the appropriate limits and reproduce the equilibrium distributions to some higher order in the integration time step. Comparisons with other existing integration schemes are carried out both for static and dynamical quantities.

Dynamic importance sampling for the escape problem in nonequilibrium systems: observation of shifts in optimal paths.

The activation problem is investigated in two-dimensional nonequilibrium systems. A numerical approach based on dynamic importance sampling (DIMS) is introduced. DIMS accelerates the simulations and allows the investigation to access noise intensities that were previously forbidden. The escape path is observed to be shifted compared to a heteroclinic trajectory calculated in the limit of zero-noise intensity. A theory to account for such shifts is presented and shown to agree with the simulations for a wide range of noise intensities.

Drastic facilitation of the onset of global chaos.

We show that the onset of global chaos in a time periodically perturbed Hamiltonian system may occur at unusually small magnitudes of perturbation if the unperturbed system possesses more than one separatrix. The relevant scenario is the combination of the overlap in the phase space between resonances of the same order and their overlap in energy with chaotic layers associated with separatrices of the unperturbed system. We develop the asymptotic theory and verify it in simulations.

Fast monte carlo simulations and singularities in the probability distributions of nonequilibrium systems.

A numerical technique is introduced that reduces exponentially the time required for Monte Carlo simulations of nonequilibrium systems. Results for the quasistationary probability distribution in two model systems are compared with the asymptotically exact theory in the limit of extremely small noise intensity. Singularities of the nonequilibrium distributions are revealed by the simulations.

Generic noise-enhanced coding in neuronal arrays.

We demonstrate that, in a parallel array of model neurons, the optimizing influence of internal noise on the global information is far greater than that reported for a single neuron. In particular, stochastic resonance (SR) effects, that optimize information transmission, occur independent of stimulus level or the setting of the neural threshold. We further show that adjusting the threshold to maximize information transmission does not remove SR effects. Consequently, and in contrast to a single neuron, in neuronal arrays noise appears to be an essential element of an optimal coding strategy.

Strong enhancement of noise-induced escape by nonadiabatic periodic driving due to transient chaos.

We have found a mechanism by which a moderately weak nonadiabatic periodic driving may significantly facilitate noise-induced interwell transitions in an underdamped multiwell system. The mechanism is associated with the onset of a homoclinic tangle in the noise-free system: if the ratio of the driving amplitude A to the damping gamma exceeds a critical value approximately 1, then the basins of attraction of the linear responses related to different wells are mixed in a complex manner in some layer associated with the separatrix of the undriven nondissipative system, and the minimal energy in such layer is lower than the top of the barrier. Thus the energy to which the system needs to be activated by the noise, to be able to make a transition, is lower than the top of the barrier.

Short time scales in the Kramers problem: a stepwise growth of the escape flux.

We prove rigorously and demonstrate in simulations that, for a potential system staying initially at the bottom of a well, the escape flux over the barrier grows on times of the order of a period of eigenoscillation in a stepwise manner, provided that friction is small or moderate. If the initial state is not at the bottom of the well, then, typically, some of the steps transform into oscillations. The stepwise/oscillatory evolution at short times appears to be a generic feature of a noise-induced flux.

Activated escape of periodically driven systems.

We discuss activated escape from a metastable state of a system driven by a time-periodic force. We show that the escape probabilities can be changed very strongly even by a comparatively weak force. In a broad parameter range, the activation energy of escape depends linearly on the force amplitude. This dependence is described by the logarithmic susceptibility, which is analyzed theoretically and through analog and digital simulations. A closed-form explicit expression for the escape rate of an overdamped Brownian particle is presented and shown to be in quantitative agreement with the simulations. We also describe experiments on a Brownian particle optically trapped in a double-well potential. A suitable periodic modulation of the optical intensity breaks the spatio-temporal symmetry of an otherwise spatially symmetric system. This has allowed us to localize a particle in one of the symmetric wells. (c) 2001 American Institute of Physics.

Noise-induced escape on time scales preceding quasistationarity: New developments in the Kramers problem.

Noise-induced escape from the metastable part of a potential is considered on time scales preceding the formation of quasiequilibrium within that part of the potential. It is shown that, counterintuitively, the escape flux may then depend exponentially strongly, and in a complicated manner, on time and friction. (c) 2001 American Institute of Physics.

Atomic micromotion and geometric forces in a triaxial magnetic trap

Nonadiabatic motion of Bose-Einstein condensates of rubidium atoms arising from the dynamical nature of a time-orbiting-potential (TOP) trap was observed experimentally. The orbital micromotion of the condensate in velocity space at the frequency of the rotating bias field of the TOP was detected by a time-of-flight method. A dependence of the equilibrium position of the atoms on the sense of rotation of the bias field was observed. We have compared our experimental findings with numerical simulations. The nonadiabatic following of the atomic spin in the trap rotating magnetic field produces geometric forces acting on the trapped atoms.

Ratchet driven by quasimonochromatic noise

The currents generated by noise-induced activation processes in a periodic potential are investigated analytically, by digital simulation and by performing analog experiments. The noise is taken to be quasimonochromatic and the potential to be a smoothed sawtooth. Two analytic approaches are studied. The first involves a perturbative expansion in inverse powers of the frequency characterizing quasimonochromatic noise and the second is a direct numerical integration of the deterministic differential equations obtained in the limit of weak noise. These results, together with the digital and analog experiments, show that the system does indeed give rise, in general, to a net transport of particles. All techniques also show that a current reversal exists for a particular value of the noise parameters.

Fluctuations and the energy-optimal control of chaos

The energy-optimal entraining of the dynamics of a periodically driven oscillator, moving it from a chaotic attractor to a coexisting stable limit cycle, is investigated via analysis of fluctuational transitions between the two states. The deterministic optimal control function is identified with the corresponding optimal fluctuational force, which is found by numerical and analog simulations.